Spin textures in slowly rotating Bose-Einstein Condensates 
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Slowly rotating spin-1 Bose-Einstein condensates are studied through a variational approach based 
upon lowest Landau level calculus. The author finds that in a gas with ferromagnetic interactions, 
such as 87 Rb, angular momentum is predominantly carried by clusters of two different types of 
skyrmion textures in the spin-vector order parameter. Conversely, in a gas with antiferromagnetic 
interactions, such as 23 Na, angular momentum is carried by 7r-disclinations in the nematic order 
parameter which arises from spin fluctuations. For experimentally relevant parameters, the cores of 
these 7r-disclinations are ferromagnetic, and can be imaged with polarized light. 
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One of the most remarkable features of recent experi- 
ments on quantum degenerate alkali atoms has been the 
observation of quantized vortices in rotating Bose con- 
densates of spin polarized 87 Rb and 23 Na [I). These vor- 
tices are a consequence of the irrotational nature of the 
superfluid flow in a scalar condensate (ie. one without 
any spin degrees of freedom). If these same experiments 
were conducted in the absence of a magnetic field, the 
velocity of these spin-1 atoms would no longer be con- 
strained to be irrotational; however, the curl of the ve- 
locity field would be fixed by spatial variations in the 
direction of the atomic spins The vortices would 
therefore be replaced by intricate spin textures. Here 
we predict the detailed spin patterns which will be found 
in such a rotating gas of spin-1 bosons. 

Our first understanding of these structures came from 
Ho [3] and Machida and Ohmi who independently 
proposed that spin textures replace vortices in spinor 
gases. Subsequently several other authors investigated 
the detailed structure of these textures in both the slowly 
rotating p| and fast rotating limits @ . The present work 
differs from these previous calculations in two vital ways: 
(i) it makes use of very simple variational wavefunctions 
from which the essential physics can be easily extracted, 
and (ii) it emphasizes the rotational properties of the 
local order parameter (which has both a vector and a 
nematic component.) The advantage of this approach, 
which explicitly considers the order parameter symmetry, 
is that it provides a scheme for classifying the spin tex- 
tures, illustrates their structure, and most importantly, 
suggests a means for experimentally detecting them. 

Although lacking our emphasis on symmetry argu- 
ments, similar theoretical techniques have recently ap- 
peared in an excellent preprint by Reijnders et al. pj. 
That preprint surveys the properties of rotating spin-1 
bosons for all possible parameters. As such, it is comple- 
mentary to the more focused study that we present. 

Throughout this paper, we work within the mean-field 
lowest Landau level approximation. This variational ap- 
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proach, detailed in section lITll is quantitatively accurate 
when the gas is so dilute that the interaction energy per 
particle is small compared to Tiu), where Ti is Planck's con- 
stant, and lj is the frequency of the harmonic trap which 
confines the gas. Our qualitative predictions, especially 
those based upon symmetry arguments, can be applied 
to denser systems. Additionally, we restrict ourselves to 
the slowly rotating limit where the angular momentum 
per particle is on the order of a few fi. In this limit, 
the angular momentum is carried by a small number of 
'elementary' textures, each of which is analogous to a 
single vortex. We reserve the discussion of faster rota- 
tion speeds, where a regular lattice of these textures are 
found, to a future Article. 

This Article is organized into the following topics: (i) a 
brief introductory example which illustrates the typical 
structure of spin textures; (ii) a review of the Hamil- 
tonian and order parameter of a spin-1 gas of neutral 
Bose atoms; (hi) a description of our main theoretical 
approach; (iv) results; and (v) broader implications. 



I. INTRODUCTION 

In this section we discuss the two component Bose gas 
in order to develop intuitive understanding of how in- 
ternal degrees of freedom affect a rotating condensate. 
Later we will apply this intuition to the spin-1 gas. In 
addition to simplicity, the two component gas has the 
advantage of extensive experimental investigations. As a 
concrete example one can consider two hyperfine states 
of 87 Rb in a magnetic trap, as produced at JILA ||. The 
order parameter here is a two-component single-particle 
wavefunction \& = (ipo > ^l) j where i/jq and ipi are the 
macroscopically occupied spatial wavefunction of each of 
the hyperfine states. Using a phase-imprinting technique, 
the experimentalists are able to place one component, 
tpi into a vortex state, while leaving ipo in the ground 
state of the trap 0. As a simple variational calculation 
shows |l0|, this unusual 'half- vortex' state is actually the 
ground state of the system for some rotation speed, and 
in the limit of weak interactions in two-dimensions can 
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FIG. 1: Top view and perspective view of a 4-7T skyrmion. At 
the center (edge) the spins point out of (into) the page, while 
half way between the center and edge the spins lie in the plane 
of the page, rotating once as one circles the center. The spins 
are located on a two dimensional plane, but point in three 
dimensions. 



be written as 
# = 



0o 
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(1) 



where z = x + iy is the coordinate in the plane, and d = 
^Jh/mui is the oscillator length, formed from Planck's 
constant H, the particle mass m and the frequency u> of 
the harmonic trap which is confining the particles. One 
can interpret \P in terms of a pseudo-spin with polar angle 
9 and azimuthal angle 0, by writing 



* = Vpe 



| cos(0/2) e"^/ 2 
sin(0/2) e^/ 2 



(2) 



where p and \ represent the local density and phase. The 
pseudo-spin behavior is then shown in figure^ In partic- 
ular, at the center of the cloud, the ipi component van- 
ishes, and the pseudo-spin points "up." At the edge of the 
cloud ipo vanishes, and the pseudo-spin points "down." In 
between, the pseudo-spin rolls from up to down, covering 
47T steradians of the sphere. This texture is often referred 
to as a skyrmion in analogy to the work of T. H. Skrymi 



[Tl| . A detailed analysis of this two component Bosc 
gas shows that these skyrmion textures are ubiquetous, 
and that in the rapidly rotating limit angular momen- 
tum is carried by a lattice of skyrmions fl2l Il3j , where 
the geometry of the lattice is sensitive to the interaction 
parameters. 

This (pseudo)spin-l/2 example teaches us that a rotat- 
ing spinor condensate tends to carry angular momentum 
by twisting its spinor order parameter. Leanhardt et al. 
at MIT [lj| recently created the spin-1 and spin-2 anal- 
ogy of this texture. Similar structures have been seen 
in rotating superfluid 3 He-A |l5j . and in quantum hall 
systems [l6| . The remainder of this Article predicts the 
equilibrium spin textures which will be found in rotating 
spin-1 condensates. 



II. SYMMETRIES OF THE ORDER 
PARAMETER 

In this section we describe the order parameter of a 
spin-1 Bose gas. 

As first discussed by Ho Q and by Machida and Ohmi 
0, a trapped gas of spin-1 bosons interacting via a 
pairwise short-range interaction which is invarient under 
global spin rotations has a Hamiltonian of the form 



H 




1 \ - 



(3) 
(4) 



where i,j running from 1 to AT label the particles, pi, 
and Si are the momentum, position, and spin operators 
for each particle, m is the particle mass, U (r) = mw 2 r 2 /2 
is the harmonic trapping potential, and V is the pairwise 
interaction which is parameterized by two constants cq 
and C2, which represent density and spin interactions. 
The most commonly studied alkali atoms have | C2 1 / 1 Cq | ~ 
0.05 <C 1. Ferromagnetic (c2 < 0) and antifcrromagnetic 
(02 > 0) interactions are respectively found in Rb and 
23 Na 0. 

In a Bose condensate of spin-1 atoms, the order param- 
eter is the three component wavefunction (ipi, 0>o> 0>_ i) 
representing the single-particle state which is macroscop- 
ically occupied. The indices 1,0,-1 represent the spin 
projection along the z direction. As pointed out by 
Machida and Ohmi 4], it is convenient to introduce 
a Cartesian representation of this wavefunction tp x = 
(-01 - V-i)/v^,^j/ = i(tpi + ip-i)/V2,tp z = -0o, so that 
?A = {Tpxjipyi Tpz) transforms as a vector under spin rota- 
tions. The order parameter is conventionally normalized 
so that the density is p(r) = ip\r)-il>[r). In the remainder 
of this paper we will freely move between the Cartesian 

[ip = (ipx,^y,ipz)] and spherical [(01,00,0-1)] represen- 
tation of the wavefunction, depending upon which is most 
useful. In general, Latin indices (such as a, b) will refer to 
the Cartesian components, while Greek indices (such as 
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fi, v) refer to the spherical components. Standard vector 
notation will be used when working within the Cartesian 
representation. 

The local expectation value of the spin in the con- 
densate is (S(r)) = (£\ 8(r - rj )Sj) = $*{r) x ^{r)/i. 
To minimize the spin interaction term for a ferromag- 
netic/antiferromagnetic gas one maximizes/minimizes 
| S\ 2 = x ip\ 2 . Consequently a ferromagnetic gas 
prefers an order parameter with the structure 



tp = n + im, 



(5) 



where n and m are real vectors with n _L m and 
|m|, corresponding to the spin with magnitude 
|^| 4 pointing out of the plane in which n and 



\n\ = 
\S\ 2 



m lie. When restricted to this maximally polarized sec- 
tor, the order parameter of the ferromagnetic gas can be 
taken to be the orthogonal triad n, m, S, whose transfor- 
mational properties are isomorphic to the group SO(3). 
Conversely, in the antiferromagnetic gas, the order pa- 
rameter prefers to have \S\ 2 = 0, which gives 







= e^n, 



(6) 



where n is a real vector and <f> is a phase. Thus antiferro- 
magnetic interactions lead to an order parameter space 
isomorphic to 52 x U(l)/Z2, where £2 is the sphere on 
which n lies, U(\) is the symmetry of the phase </>, and 
the quotient with the two-element discrete group Z2 rep- 
resents the fact that ip is invariant under simultaneously 
reversing the direction of n and taking cf> — > <f) + ir. With- 
out the phase, the order parameter space is the group 
S2/Z2, which can be identified with the projective plane 
RP2- The importance of this projective structure has 
been emphasized by Zhou |l7j . 

A more intuitive approach to understanding the sym- 
metries of the spin-1 condensate comes from considering 
the moments of the local spin operator. A ferromagnetic 
interaction favors a maximally polarized state; for exam- 
ple (■01, ipo, ip-i) = (1,0,0), for which the local spin is 
(S) = z. In this case one can take (S), an object which 
transforms as a vector, to be the relevant order param- 
eter for discussing spin textures. Antiferromagnetic in- 
teractions favor a minimally polarized state; for example 
(V>i, -00,-0-1) = (0,1,0), for which (S) =0. To find a use- 
ful order parameter, one must consider the fluctuations 
of the spin 5S a SS b = Re{S a S b ) - (S a )(S b ) [3, which 
for the (0, 1,0) state is SS a SSb = S a b(l — S az ), where 5 a b 
is the Kroneker delta. Consequently, this state has no 
spin fluctuations in the z directions, and isotropic fluctu- 
ations in the x-y plane. Due to the azimuthal symmetry, 
the spin fluctuation tensor 8S a 8Sb maps onto itself un- 
der rotation by 7r radians. Equivalently, one can say that 
the spin fluctuations represent a nematic order parameter 
(meaning that it transforms under rotation as an arrow 
with no head on it.) 

As a consequence of the disparate types of spin tex- 
tures which involve vector and nematic order parameters, 
gases with ferromagnetic and antiferromagnetic interac- 
tions behave in quite distinct manners. 



A. Local Vorticity 

We now discuss how the symmetries of the order pa- 
rameter are related to the curl of the velocity field, 



2im 



(7) 



where repeated indices are summed over. The key obser- 
vation is that unlike a scalar condensate, the spin-1 con- 
densate can have a velocity field with extended vorticity. 
This property is most simply understood by introducing 
a tensor 



(8) 



which carries all information about the local spin or- 
der parameters but no information about the density or 
phase. By construction Q is Hermitian (Q a b — Qta) ano - 
has trace 1. We decompose Q into a sum of irreducible 
tensor operators, 



Qab 



Q {0) 5 ab /3 + ie abc Q 
1, 

s=(S)/p, 
(2/3)5 ab - ((S a S b ) 



(i) 



/2- 



)(2) 
tab ' 



(S b S a ))/2p 2 



(9) 
(10) 

(11) 
(12) 



where e a bc is the totally antisymmetric unit tensor and 
p is the density. The three irreducible components are 
a constant scalar Q^°\ a vector Q^ 1 ) which represents 
the local spin order, and a symmetric traceless tensor 
Q( 2 ) which represents the local spin fluctuations. The 
magnitude of these various components are constrained 
by the fact that Q is idempotent (Q 2 — Q), and hence 







= Tr(Q 2 ) 
= -2/3-1 



- (TrQ) 

Q( 1 )-Q( 1 )/2 + Tr((Q( 2 )) 2 ). 



(13) 
(14) 



The last two terms are non-negative scalars which rep- 
resent the amount of vector and nematic order. This 
expression clearly shows that these are competing order 
parameters, as increasing one of them requires reducing 
the other. 

As demonstrated in appendix ^ the vorticity can be 
written in terms of Q as 

Ti 

V X V S = i — Qabi^Qbc) x (VQ C „). (15) 

m 

In the spin polarized case 2Q a b = Sn b —s n Sb+icnhr.Sr. (with 
s • s = 1) one recovers the Mcrmin-Ho [2(J relationship 



V x v s 



-tabcSa(VSb) X (Vs c ). 



(16) 



In the purely (uniaxial) nematic case Q a b = n a rib 
(with n - n = 1), one instead finds an irrotational flow, 
V x v s = 0. Note that equation (|15fl is correct for all 
spin, not just spin 1. 
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In a more geometric language, the superfluid velocity 
is the phase one-form associated with ip, its curl is the 
phase two-form, and Q is the operator which projects 
onto ip- The relationship (f 1 51) between the projector and 
the phase two-form, is generally valid, and shows up in 
other contexts [Tflj . 

B. Visualization 

The ferromagnetic gas has a vector order parameter 
(S) which is visualized by drawing a vector at each point 
in space (as in figure^. The antifcrromagnetic gas has 
a more complicated order parameter. In the strongly 
antiferromagnetic limit, the order parameter has the form 
of eq. J?SJ|, and one can visualize the order by plotting a 
rod aligned along the direction n at each point in space. 
In general, the order parameter deviates from eq. iJBJ, 
and it is more convenient to describe the order in terms 
of Q< 2 ) or QW = (Q afc + Q ba) / 2 = (1/3)<U + Q$. The 

tensor is simpler to work with than as it is 
positive semidehnite (meaning that all eigenvalues are 
greater or equal to zero). In the case where ip is given 
by then has a single non-zero eigenvalue whose 
eigenvector points in the n direction. In the general case, 
one can define h to coincide with the direction of the 
largest eigenvector of Q^ s \ 

Generically has three distinct eigenvalues, and 

therefore describes a biaxial nematic, which cannot be 
represented solely in terms of n. For a more complete 
picture of the local nematic order, one replaces the rods 
by elipsoids whose three principle axes coincide with the 
eigenvectors of Q( s >. The length of the principle axes are 
taken to be proportional the eigenvalues. These ellip- 
soids are easily constructed by treating as a linear 
transformation to a sphere. Starting from a sphere S, 
the resulting ellipsoid is defined as the set of points sat- 
isfying r' a — Q^ r bi f° r r S S. In the uniaxial limit, 

(s) 

Qab = Qab = n a nb, two of the ellipsoid's axes have 
length zero, and this degenerate ellipsoid becomes a rod 
pointing in the direction n. In the ferromagnetic limit 

( s) 
2 Qab 

— ®ab Sa Sb, the ellipsoid becomes a disc, whose 
normal is the spin vector s. 

In both the ferromagnetic and antiferromagnetic gas, 
the spins and directors point in three dimensional space. 
Due to this intrinsic three dimensionality, it is impos- 
sible to faithfully represent the spin textures on a two- 
dimensional sheet of paper. Consequently, in this Article 
we only provide sketches of the simplest textures. Ani- 
mated representations of the more complicated structures 
can be found in the EPAPS archive |2lj . 

For simplicity, wherever possible we attempt to de- 
scribe the nematic textures solely in terms of n. 

III. LOWEST LANDAU LEVEL CALCULUS 

In this section we present a calculational scheme for 
studying Eq. ©. Our approximations are motivated 



by the limit where interactions are weak compared to 
the trapping potential. Weak interactions are naturally 
reached under fast rotation, as in the rotating frame the 
centrifugal force effectively reduces the trap strength, 
causing the cloud to spread out over a larger area, reduc- 
ing the density and hence the interactions. Although mo- 
tivated by this ultra-dilute limit, our approach gives an 
effective approximation to the properties of a system with 
much stronger interactions. For example, both static [2^ 
and dynamic |23T | properties of vortex lattices in a scalar 
condensate have been very successfully modeled by this 
method. 

In the rotating frame the Hamiltonian is shifted to 
H' = H — fi ■ L, where the vector fl points along the axis 
of rotation with modulus given by the rate of rotation, 
and L is the angular momentum vector of the system. 
We restrict ourselves to two dimensions (taken to be the 
x-y plane), and rotate about z axis. This dimensional re- 
striction is consistent with the weakly interacting limit, 
where only a single mode of the trapping potential is 
occupied in the £ direction. To extend this method to 
describe the more general (3D) situation, one can use a 
Thomas- Fermi approximation as in |'2l|. 

In the absence of interactions, the eigenstates of 
the two-dimensional single particle Hamiltonian H' sp = 
p 2 /2m + mw 2 r 2 /2- Q,(r x p y — r y p x ) have energies E nm = 
hcu + h(u> + f2)n + h(u> — fi)m, and angular momentum 
Lnm = h(n — rn), with 7i, m = 0, 1, . . . If the interac- 
tion energy per particle is sufficiently small compared to 
h{ui + CI) then one can build up the quantum state from 
the "Lowest Landau Level" (LLL) , made up of the states 
with n = 0. Here we restrict our analysis to the mean 
field level, where one single particle state is macroscopi- 
cally occupied. The most general single particle state in 
the LLL can be written as 



^{x,y)=Y,^^=e-^' 2 , (17) 
\Jttv. 

3 



where z = (x + iy) /d is the coordinate in the plane in 
complex notation, scaled by the oscillator length d = 
y/h/mui. The coefficients Cj^ can also be expressed in 
Cartesian notation Cj — (cj x , Cj V , Cj Z ). We normalize this 
wavefunction so that 



Jdxdy ■tf=J2 c *r c i = 1 - ( 18 ) 



Substituting this ansatz into @ , and neglecting the zero 
point energy, the scaled energy per particle in the rotat- 
ing frame becomes 
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(H - QL Z ) 
Nh(uj - fi) 

1 

Ti{w - fi) 

E [c? + m 



(19) 



2m 



2 'XKiV 2| '" , - ( '- ! " 



[(l + c 2 )(Trd;.d j )-c 2 |Trd j | 2 ] 



2tt(w - fi)' 

r - 



Here 'Tr' signifies taking the trace of a matrix. The 
simplicity of these equations illustrate the utility of the 
Cartesian representation of ip. In the LLL, the system is 
parameterized by two dimcnsionlcss couplings; r), which 
measures the ratio between the interaction strength and 
the level spacing, and c 2 = c 2 /co, the relative strength 
of spin and density interactions. Larger r\ corresponds to 
faster rotation. If one rotates at speeds where Q > u), 
the centrifugal force becomes stronger than the trapping 
force and the system becomes unstable. This centrifu- 
gal limit coincides with rj — ► oo. One can experimentally 
tune 77 between and +00. 

In this Article we find the ground state as a func- 
tion of 77 and c 2 by minimizing l|20l) with respect to the 
parameters cj using a conjugate gradient method. We 
truncate the allowed values of j to be j = 1, 2, • • • , J. 
Convergence is ensured by taking several values of J, 
ranging from 2 to 128, and finding that for sufficiently 
large J the energy saturates. Although this variational 
approach is much simpler than directly minimizing the 
Gross-Pitaevskii equation, it is still a nontrivial numeri- 
cal task. Since cj is a three component complex object, 
one has the equivalent of 6 J real variational parameters. 
Symmetries can be used to reduce the number of inde- 
pendent parameters. All our results use the approximate 
physical value, c 2 = ±0.05 [2{|. For a discussion of other 
values of c 2 , see @- 

Not surprisingly, the energy landscape is quite com- 
plicated with many metastable local minima. Computa- 
tionally, this large degree of metastability makes finding 
the absolute ground state quite difficult. In section lTVI we 
show how symmetries can be used to help find these min- 
ima. This symmetry approach is augmented by a brute- 
force search strategy where we run our minimization rou- 
tines from a very large number of initial conditions, and 
select out the lowest energy minima found. Experimen- 
tally, the presence of so many local minima implies that 
the state found in the laboratory will depend strongly 
upon the method of preparation. Moreover, one would 
expect to see domain structures where different parts of 
the sample exhibit different phases. 



A. A More General Ansatz 

Although all the numerical calculations presented here 
involve equation (|20|) . it should be noted that by intro- 
ducing an extra variational parameter, one can greatly 
extend the realm of validity of the LLL approximation. 
The resulting equations are not much more complicated 
than <HIJ|. 

The extended LLL ansatz scales the lengths in (|17f) by 
a new variational parameter v so that z = v 1 / 2 (x + Iy)/d. 
Neglecting zero point motion, the energy per particle is 



Nh 

h 

I? 







- n 













h + (u-Sl)juI 2 , (22) 



E(i + 1 ) c r c J. ( 23 ) 

j 

^[(l + caXTrdr.djJ-^lTrdjl 2 ]. (24) 



Minimizing with respect to v with {cj} fixed yields 



H' = Jift 



U)\ 1+1]- 



1-2 U) — 



- n 



1 + Vh u-Sl 
2/i to ' 



(25) 
(26) 



In the limit A = r)[(oj — ri)/w][/ 2 //i] -C 1 one recovers 
(|20|l . One can minimize l|25l) by the methods already 
discussed. Note that instead of depending upon only 
two experimental inputs, 77, and c 2 , the minima in the 
extended LLL ansaztz also depends on the dimensionless 
ratio (lu — f2) /ui. 



IV. SYMMETRIES OF THE SPIN TEXTURES 

A. Symmetries 

Given the complexity of the energy landscape, an unre- 
stricted minimization of eq. I|20|) is a daunting task. The 
large number of variational parameters makes the com- 
putation expensive, and the large number of metastable 
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FIG. 2: Top view and perspective view of an 87r skyrmion. 
This texture differs from figureQin that the planar projection 
of the spins rotates twice as one circles the origin. 

minima means that an exhaustive search is needed to 
find the absolute minimum. Moreover, a classification 
scheme is needed to describe the minimal configurations 
once they are found. We solve all of these problems by 
studying the possible symmetries of the spin textures. In 
section IVIES! we discuss how these symmetries can be 
used to detect the textures. This strategy was inspired 
by discussions with Dan Rokhsar, who, with Dan Butts, 
used similar techniques in understanding vortex lattices 
in scalar condensates p2|. 

Our Hamiltonian @ is separately invarient under each 
of the following operations: (i) rotating spatial coordi- 
nates about the origin by an arbitrary angle 9, 

4>{r a ) — ► R ab (9)^(r b ) (27) 

R(0) = ( C ° S l Sin ^; (28) 
w V — sin 6 cos 1' v ' 

(ii) Rotating all spins by arbitrary Euler angles <fi, 9, and 
X, 

Vv — > ?WM,xM, (29) 

K{9,dp, x ) = e-txS'e-HS'e-** 8 '] (30) 

where S X ,S V ,S Z , are the operators for the spin com- 
ponents; (iii) simultaneously reflecting spatial coordi- 



nates across the line perpendicular to the unit vector 
n, and applying the time reversal operator (which takes 

V^(r) — > ^%(R fl r) (31) 
R fl r = r-2nr; (32) 

(iv) reflecting the spin across the plane perpendicular to 
the unit vector n, 

— ► Knlf (33) 
Hn^f = ^f-2h-^. (34) 

and (v) a global gauge transformation by an arbitrary 
phase x, 

VV — » e ix V M - ( 35 ) 

Although the Hamiltonian is invarient under each of 
these transformations, the mean-field wavefunction is 
not. We will classify spin textures by the way in which 
they break these symmetries. 

B. Scalar 

To see how these symmetries manifest in spin textures, 
it is helpful to first consider a scalar condensate [2^, for 
which only symmetries (i), (hi), and (v) are relevant. In 
the lowest Landau level approximation, the non-rotating 
scalar ground state is ipir) = e~' r ' , which is invarient 
under (i) and (iii), but not (v): i.e. the only broken 
symmetry in the non-rotating condensate is the gauge 
symmetry. 

At appropriate rotation speeds the ground state con- 
tains a single vortex, and (in the lowest Landau level) the 
order parameter becomes ip = ze~' r ' / 2 , where, as usual, 
z = x + iy. This single vortex state is no longer invarient 
under any of the transformations, (i), (iii), or (v). How- 
ever, it is invarient under appropriate combinations. For 
example, rotating the spatial coordinate by any angle 9 
transforms z — > e lS z, and can therefore be undone by a 
global gauge transformation l|35|) with \ = 9. This is an 
example of a continuous symmetry, as for each 9 there 
is a combined spatial rotation and gauge transformation 
which leaves the state invarient. This symmetry is not re- 
stricted to the lowest Landau level approximation, and it 
should be clear that any condensate with a single vortex 
at the center will have this symmetry, expressed as 

</>(R(0)r) = e* OT V(r), (36) 

where m is a fixed integer which gives the number of 
quanta of circulation which are concentrated in the vor- 
tex core. An ordinary vortex has m = ±1. 

At faster rotation speeds, the condensate will contain a 
small cluster of vortices. As an example, one can imagine 
two vortices, symmetrically placed at z = T^Oi so that 
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the order parameter is -0 = (z — z )(z + z )e~l r l 2 / 2 . No 
continuous symmetry exists here. However, there does 
exist a discrete symmetry in that this state is invarient 
under rotating space by 9 — ir. Similarly, three vortices 
which form an equilateral triangle will be invarient under 
a rotation by 9 = 2tt/3. Such discrete symmetries are 
generally of the form 

^(R(27r/a)r) = e 2 ™"/ a V(r), (37) 

where a is an integer which, together with the integer m, 
describes the symmetry of the cluster of vortices. 

The wave function for a cluster of vortices with a re- 
flection plane will be invarient under a combination of 
a reflection across a fixed line, time reversal (complex 
conjugation), and a gauge transformation, 

V>(fW)* - ±V(r) (38) 

where the line of reflection is perpendicular to the unit 
vector n. We have written the gauge transformation as 
a multiplication by ±1, as the only consistent phases are 
X = 71", 0. 

C. Spin 1 

As with vortices in a scalar gas, spin textures in a 
spinor condensate are characterized by their properites 
under different combinations of symmetry operations. 
The spin degrees of freedom only introduce additional 
symmetries. 

At low rotation speeds one expects to find a continu- 
ous symmetry, generically involving the simultateous ro- 
tation of space, spin, and phase. Without any loss of 
generality, we can take the spins to rotate about the z 
axis, in which case the symmetry is formally described by 
the statement that there exists constants m and n such 
that for all 9, 

^(R(0)r) = e" ne Tl^(n9)Mr), (39) 

where R(#) represents a spatial rotation about the origin 
by an angle 9, TZ^ V (4>) represents a spin rotation about 
the z axis by an angle 0, and m, n count the number of 
times the phase and spin angles rotate in comparison to 
the spatial angle 9. Single valuedness of the left hand 
side of l|39|) constrains m and n to either both be inte- 
gers, or both be half- integers. Further explanation of this 
constraint will be given in section llV C II 

Textures obeying (|39[) will be described as the 'elemen- 
tal,' 'single,' or 'azimuthally symmetric' spin textures. 
They are the building blocks for all more complicated 
strucures. 

At higher rotation speeds, this continuous symmetry 
is replaced with a discrete one. In analogy to the scalar 
case <|37[) . the discrete rotational symmetry is of the form 

^(R(27r/o)r) = e" n2 ^ a 1l^(n2ir/a)Mr), (40) 



where a is a fixed integer describing the geometry of the 
texture. 

In addition to rotational symmetries, the state may 
have reflection symmetries. The most general form of 
reflection symmetry which we consider consists of simul- 
taneously performing transformations (iii), (iv), and (v), 
namely reflecting the spatial coordinates accross a line, 
time reversal, reflecting the spins accross a plane, and 
multiplying by the factor ±1. By globally rotating our 
spins, we can always assume that the spin and spatial re- 
flection planes coincide. Mathematically this symmetry 
can be expressed as 

(^(Riir))* - ±(Ofi) lia ,Mr), ( 41 ) 

where h is the normal to the reflection plane, R„ is the 
spatial reflection and (TZn)^ is the spin reflection. One 
can also consider the case where no spin reflection is per- 
formed 

(^(R fl r))* = ±^W- (42) 

We will refer to (|41|l as a reflection and 142|) as a partial 
reflection. 



1. implications 

Enforcing the various possible symmetries describe in 
equations (|39|1 through l|42|l . greatly reduces the num- 
ber of parameters which must be minimized. Given the 
explict form of the rotation operators 

R{9)z = e- l9 z (43) 
n^(9) = 8^ u e ie ^ (44) 

and the lowest Landau level ansatz ip^ = 
J2j Cpj-z-'e - ' 2 ' 2 / 2 , the continuous symmetry J35J|, 
restricts c w to be zero unless 

j = m + vn. (45) 

Since v — 1, 0, — 1, this requires any non-trivial texture to 
have m, n both be integers, or both be half integers. At 
most three c M: ,'s are non-zero. By combining an overall 
spin rotation, a spatial rotation, and a gauge transfor- 
mation, these three parameters can be taken to be real. 
Using the normalization constraint X^'u \ c n\ ~ 1> one 
is left with only two (real) variational parameters. 

The discrete symmetries similarly reduce the number 
of parameters. For example the discrete rotation sym- 
metry forces c^j — unless 

j = m + vn (mod a). (46) 

Up to an overall phase, enforcing the reflection symmetry 
(EfTl) with n — y, requires c^j is real. 
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V. RESULTS 
A. Azimuthally Symmetric Textures 

In the limit of very slow rotation, the angular momen- 
tum will be carried by a single 'elemental' spin texture, 
characterized by the continuous symmetry in equation 
(|39|1 . This structure is analogous to a single vortex. At 
faster rotation speed a cluster of these textures will form. 
Eventually the clusters will become regular lattices. In 
this section the structure of the 'elemental' spin textures 
will be explored. 

1. Ferromagnetic 

The ferromagnetic gas displays two different elemen- 
tal textures, both of which are similar to the skyrmion 
encountered in the spin-1/2 case. These two textures, il- 
lustrated in figures ^and 13 differ by the amount of solid 
angle (47r or 8n) traced out by the spin vectors as they 
move from pointing straight up at the center to straight 
down on the edge, and will therefore be referred to as 4-7r 
or 87r skyrmions. If one looks down from above, as shown 
in the top panels of figures ^ and |3 one can distinguish 
these textures by the number of times which the planar 
projection of the spin rotates when one circles the origin 
(once/twice for the 47r/87r texture). Higher order tex- 
tures (where angles greater than 87r are subtended) are 
never found to have continuous azimuthal symmetries, 
and are always more appropriately described as a collec- 
tion of 47r and 87r skyrmions. At faster rotation speeds, 
angular momentum is carried by a lattice of 47r or 87r 
skyrmions where their quantization axes takes on several 
different orientations. These textures only cover a full An 
and 87r steradians of solid angle when they exist in iso- 
lation in an infinite system. In a lattice, the skyrmions 
often touch before they can completely cover the sphere. 

The 47r and 87r textures are described by l)39|) . with 
m = n = 1 and m = n = 2. In a basis with spin 
projections (1,0,-1), the spinors take the form 

\n = m = l) = \ az | e^ 2 ' 2 , (47) 



\n = m = 2) = az 2 e^' 2 . (48) 

W) 

The real coefficients a and b are given by minimizing (|20[1 . 

In the 47r skyrmion, the velocity field has its great- 
est curl at the center of the texture, while for the 87r 
skyrmion the curl is greatest on a ring of finite radius 
(see figure ISJ) . According to the Mermin-Ho relationship 
(|15fl one can attribute these two distinct behaviors to the 
fact that the spins 'bend' fastest in these regions. 

The 8n skyrmion can alternatively be interpreted as a 
composite of four An textures j2|| [27J ■ This interpreta- 




FIG. 3: Local vorticity z ■ V X v, for the An (a) and 87r (b) 
skyrmion. Darker colors represent larger vorticity. 

tion is illustrated in figure where the spins in figure [21 
are rotated around the y-axis by 90 degrees. There are 
four points where the spins are pointing into/out of the 
page. These points lie on the outer edge of the ring of 
maximal V x v, and can be taken to be skyrmion cores. 
Note that their positions are not unique, and by glob- 
ally rotating all of the spins, these points move around 
in a ring. This lack of uniqueness is a consequence of the 
high degree of symmetry of this state. The 87r steradians 
of solid angle subtended by the texture can be divided 
evenly among these smaller skyrmions, allowing us to 
attribute In steradians to each of them. Two of the sub- 
textures, labeled (A) and (C), have cores pointing out 
of the page, while the other two, labeled (B) and (D) 
point into the page. Equidistance between these cores 
the spins lie in the x-y plane. When circling (A) or (C) 
these planar spins rotate in a positive sense while around 
(B) or (D) they rotate in a negative sense. 

The An skyrmion has lower energy than the uniformly 
polarized state (i/ji,ipo,ip-i) — (1,0,0) if and only if 
77 > 2/(l + C2). This is a second order (continuous) phase 
transition. Between the region of stability of the 47r and 
8n skyrmions, a more complicated texture is found. De- 
tails of this intervening state will be given in section lVBI 
A graphical comparison of the energies of these states as 
a function of 77 is made in figure [7] 



2. Antiferromagnetic 

The antiferromagnetic gas, with its nematic order 
parameter, supports a much different set of textures 
than the ferromagnetic gas. In particular, one finds 
that angular momentum is predominantly carried by n- 
disclinations, which as illustrated in figure are objects 
around which the nematic director rotates by 180 de- 
grees {n radians). There is no consistent way to de- 
fine the direction of the director at the center of the n- 
disclination, and one therefore calls these textures "topo- 
logical" meaning that there is a loss of continuity at the 
core. The topological nature of this excitation makes it 
very similar to a vortex, for which there is no way to 
define the phase of the order parameter at the center. 
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FIG. 4: (color online) Top view and perspective view of an 87r 
skyrmion. This is the same texture as fig. El except all spins 
have been rotated by 90 degrees about the y axis. This texture 
can be interpreted as a composite of four 4ir skyrmions, whose 
cores have been marked with circles, and labeled by letters 
(A) through (D). The spins in (A) and (C) point out of the 
page and the planar projections of the spins wrap the equator 
of the order-parameter sphere in a positive sense, while (B) 
and (D) point into the page and the projected spins wrap the 
equator in a negative sense. 



FIG. 5: Representations of a tt disclination in a nematic. 
Lines represent the local orientation of the nematic director. 
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FIG. 6: Structure of the nematic ring texture (K). Horizontal 
axes represent distance from the center of trap along the x 
axis, measured in units of the trap length. Top panel shows 
the density p. At each position in the middle panel, an arrow 
is drawn which represents the direction and strength of the 
local spin. At each position in the bottom panel, a rod points 
in the direction of the nematic director n, corresponding to 
the largest eigenvalue of Q^ s h The length of these directors 
are scaled by the total density, so that their lengths are re- 
lated to the amount of local nematic order. To construct the 
full three-dimensional spin texture, one rotates this picture 
around the origin, so that the spins near the center look like 
the skyrmion in figure Q and the nematic order away from 
the center forms a crown texture. 



In the case of a vortex, this lack of continuity causes the 
density to vanish at the core. For the experimentally rele- 
vant spin-1 gases, the density interaction is much stronger 
than the spin interaction, and it is favorable to fill the 
core with ferromagneticaly ordered atoms. 

Thus, angular momentum is carried by tt disclinations 
with ferromagnetic cores. The spins in the core align 
perpendicular to the plane in which the nematic directors 
lie. In the fast rotating limit, one finds a square lattice 
of these tt disclinations with their cores aligned in an 
antiferromagnetic checkerboard pattern. 

The tt disclination is described by (|39l) with n = m = 
1/2. In a basis with spin projections (1,0,-1), the spinors 
take the form 



\n = m = 1/2) = 




where the real number b is given by minimizing H2()(l 
The 7r-disclination has a lower energy than the uniform 
nematic state (i/)x, ipoi ty-i) — (0,1,0) if and only if 
77 > 2 — 2 v / 2c2 + 0(c 2 ). This is a first order (discon- 
tinuous) phase transition. 

We find one other azimuthally symmetric texture in 
the gas with antiferromagnetic interactions. This state 
has the same mathematical structure as Q470 , except that 
here b > 0, while in the ferromagnetic case b < 0. As 
sketched in figure this texture consists of a nearly uni- 
form nematic ring with a ferromagnetic core. The spins 
in the ferromagnetic core bend like the skyrmion in fig- 
ure H The nematic ring has its directors canted slightly 
from the z axis, forming a crown-shaped texture. We 
do not find any composite textures which include this 
structure. 
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FIG. 7: (Wide) Scaled energies £ as a function of rotation speed, parameterized by rj defined in Eq. 121H . Larger 77 corresponds 
to faster rotation. Top/bottom panels show ferromagnetic/antiferromagnetic interactions with C2 = =F0.05co. Each curve 
represents a state of different symmetry, as described in the text and figures 1101 through 1111 More detailed images of these 
states can be found in the EPAPS archive |2l| Energy scalings (Top: £ = £ 2 - 1.7125??, Bottom: £ = £ 2 - 1.80847?) are chosen 
to aid in comparing these different curves. 



B. Composite Textures 

At higher rotation speeds, angular momentum is not 
carried by single 'elemental' textures, but rather by a 
small collection of these objects. At very fast rotation 
speeds one expects to find a regular lattice. We study 
these more complicated objects by minimizing (|2TJ|) . se- 
quentially constraining the wavefunction to have the 
symmetries in (I39|l through l|42|) . 

Our numerical results are summarized in figures [7| 
through ^3 Figure shows the energy of states with 
various symmetry properties. The data is scaled so as to 
bring out the important features. There are two parts 
to the scaling: (1) £ 2 is plotted rather than just £ be- 
cause in the fast rotating limit one expects £ cx 77 1 / 2 . 
This dependence is made apparent by noting how vari- 
ous terms in the energy scale with the total size of the 
cloud; I2 oc (r 2 ) and £4 oc l/(r 2 ). Minimizing with re- 
spect to (r 2 ) then gives £ oc rj 1 / 2 . (2) We subtract /3?? 
from £ 2 , where (3 is an emperically determined constant. 
Removing this linear term makes the differences in the 
curves easier to see. We use (3 = —1.7125 (—1.8084) for 
the ferromagnetic (antiferromagnetic) data. 

The curves in figure are labeled by the letters (A) 
through (S) and can be described as follows: 
Ferromagnetic States: (A) The uniformly polarized 
state. (B) The 47r-skyrmion state \n — l,m = 1) de- 
scribed in section IV A II (C) A texture with a single 
reflection plane consisting of two elementary structures: 
a 47T skyrmion near the origin, and a second object at 
whose center the ferromagnetic order vanishes, to be re- 
place by nematic order. (D) The 87r-skyrmion state 
\n = 2, to = 2) described in section IV A II (E) Simi- 
lar to (C), except the 4-7T skyrmion is replace by an 87T 
one. (F) A texture with two reflection symmetries and 
a discrete rotation symmetry (with a=2,n=m=l). Two 
87T skyrmions with their axes canted with respect to one- 
another lie on one axis, while two weakly nematic re- 
gions, where the ferromagnetic order is reduced, lie on 




FIG. 8: Angular momentum per particle L and spin polar- 
ization per particle a = |J'd 2 rS(r)| of ferromagnetic spin 
textures as a function of rotation speed, parameterized by 77. 
Both are measured in units of h. Note the discontinuities 
when the ground state symmetry changes. 
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FIG. 9: Angular momentum per particle L and spin polariza- 
tion per particle a of nematic spin textures as a function of 
rotation speed, parameterized by 77. 



the other. (G) A texture with two reflection symmetries 
and a discrete rotation symmetry (with a=4,n=l,m=0). 
At the center lies a nematic region. Four 87r skyrmions 
with canted axes form a square surrounding the center. 
(C) as with (G), but the square is distorted into a rhom- 
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bus, resulting in a lower symmetry (a=2,ri=l,m=l). (H) 
A texture with four reflection axes and a discrete ro- 
tation symmetry (a=4,n=l,m=l). There are four 87r 
skyrmions, symmetrically situated with their axes nearly 
lying in a plane. 

Nematic States: (I) The uniform nematic state. (J) 
The 7r-disclation with ferromagnetic core described in 
section IV A 21 (K) A nematic ring with a ferromagnetic 
core described in section IV A 21 (L) Four 7r-disclinations 
forming a square. The feromagnetic cores are aligned an- 
tiferromagnetically. This state has two reflection planes, 
and a four-fold rotational symmetry (a=4,n=2,m=2). 
(I/) As with (L), but the square is slightly deformed into 
a rhombus, and the perfect antiferromagnetic alignment 
of the cores is slightly distorted. The rotational sym- 
metry is reduced (a=2,n=l,m=l). (M) Five tt disclina- 
tions form a regular pentagon. The ferromagnetic order 
at the cores lies completely in the x-y plane. Contains 
a reflection plane and a five-fold rotational symmetry 
(a=5,n=2,m=0). (M') As with (M), but the perfect five- 
fold symmetry is slightly distorted with the pentagon of 
disclinations stretched along one axis. This state only 
has a single reflection plane and no rotational symme- 
try. (N) Dominated by two stripes of three tt disclina- 
tions organized into a (nearly) square lattice, this tex- 
ture has a full reflection plane and a partial reflection 
plane (accross which the spatial but not spin degrees of 
freedom arc reflected.) (O) Eight tt disclinations: one 
at the center, surrounded by seven others. The cen- 
tral core points in the z direction, while the surround- 
ing cores are canted slightly in the — z direction from 
the x — y plane. Contains one reflection plane and a 
seven-fold rotational symmetry (a=7,n=3,m=-3). (P) 
Nine it disclinations forming a square lattice with their 
ferromagnetic cores aligned antiferromagnetically. This 
state contains two refection planes and a four-fold rota- 
tion symmetry (a=4,n=l/2,m=l/2). (P') As with (P), 
but the square lattice is deformed towards having a tear- 
drop shaped set of eight disclinations surrounding the 
central one. This state only has a single reflection plane 
and no rotational symmetries. (Q) A pattern of ten 
disclinations; two central disclinations surrounded by a 
distorted oval of eight others. This state has a single 
reflection plane. (R) A distorted square lattice, consist- 
ing of twelve disclinations. Contains one reflection plane. 
(S) A distorted square lattice of fourteen disclinations. 
Contains two reflection planes, and a discrete rotation 
symmetry (a=2,n=2,m=2). 

Three dimensional animated representations of these 
states are stored on the epaps archive [2l| . 

From this data one can calculate several observables 
including: the angular momentum and the degree of spin 
polarization as a function of r\. These results are shown 
in figures |S] and [5J The angular momentum is measurable 
through collective mode experiments |28| , while spin po- 
larization could be measured through magnetic suscept- 
ability. In figurcs^|andEl we show spatial distributions 
of the density p, vorticity z • V x v, spin density \S\ 2 , and 
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FIG. 10: Spatial structure of spin textures in a gas with 
ferromagnetic interactions. From left to right, the columns 
represent density p, vorticity z • V x v, spin density |S| 2 , 

and nematic order Q^bQ^a — P^Q^bQ^ba ■ Darker colors 
corresponds to larger magnitudes. Although not universally 
true, 4-7T (87r)-skyrmions tend to show up as black (white) dots 
in column 2 and white (black) dots in column 3. 



nematic order Q ab Q ba for each of these textures. 



VI. DISCUSSION 

A. Alternative classification of skyrmions 

The division of a given texture into elemental 
skyrmions is not unique, as is illustrated by the two 
ways of thinking about the 8tt skyrmion in sect ion IV A II 
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FIG. 11: Spatial structure of spin textures in a gas with 
antiferromagnetic interactions. With the exception of texture 
K, all black dots in the third column (which also coincide with 
white dots in the fourth column) can be identified as the cores 
of 7r-disclinations. 



Above we chose to discuss composite textures in terms of 
the spin behavior at the density maxima (which is anal- 
ogous to looking at the 87r skyrmion from the angle in 
figure |2l ■ This is a natural taxonomy scheme in that it 
reduces the total number of skyrmions which need to be 
considered, at the cost of introducing 8ir skyrmions. It 
is enlightening to also describe the texture in terms of 
the spin behavior at the density minima (analogous to 
the view of the 8n skyrmion in figure 0] in terms of four 
smaller textures). The advantage of this latter viewpoint 
is illustrated by considering texture (H). At the local min- 
ima of the density, the spin points into or out of the page 
in a manner suggestive of two interpenetrating square 
lattices. The textures around each minima can be de- 
scribed as Mermin-Ho skyrmions (or merons) pcj where 
at the center the spin points in the ±z direction, and 
then rolls over to lie in the plane. In the vicinity of each 
texture the spins trace out 27r steradians, covering half 
of a sphere. 

Under this interpretation the following states have new 
descriptions: (D) Four skyrmions form a square. Two of 
their quantization axes point into the page, two out of 
the page. (E) Along the x-axis the density has three 
minima, and there are correspondingly three elemental 
textures, two skyrmions and a nematic region. (F) Three 
skyrmions lie along the x-axis. (G) Four skyrmions form 
a square and a nematic region sits at the center. (C) 
As (G) but the skyrmions form a rectangle. (H) Two 
interpenetrating square skyrmion lattices with four 4ir 
skyrmions pointing into the page and five pointing out. 



B. Experimental Consequences 

In this section we address the questions of how to ex- 
perimentally create and measure these spin textures. 

Creation: One should be able to create these spin 
textures by the same techniques used to create vortex 
lattices in scalar condensates yj. These methods include; 
'stirring' the cloud with a detuned laser, rotating an el- 
lipsoidal trap, and cooling a rotating cloud through the 
BEC phase transition. Several caveats must be kept in 
mind however: (i) stray magnetic fields and magnetic 
field gradients must be minimized, and (ii) the large de- 
gree if degeneracy in this system combined with exper- 
imental randomness may lead to more complicated spin 
textures than those seen here. In particular, at fast rota- 
tion speeds one would expect to find a domain structure, 
where different parts of the clouds contain different lat- 
tices of skyrmions/disclinations. 

Detection: We propose directly imaging these tex- 
tures. First, one turns off the atomic trap, allowing the 
cloud of atoms to expand. If a magnetic field gradient 
is introduced during the expansion, the different compo- 
nents (ipi,ipOiip-i) w iU become spatially separated as in 
a "Stern-Gerlach" experiment. Each of the three com- 
ponents can then be separately imaged. This detection 
method has been used to observe skyrmion textures in 
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FIG. 12: Densities ?/<i| 2 , | Vo | 2 , and IV'-il 2 of the three compo- 
nents of texture (H), as seen from two different quantization 
axes. Darker colors represent higher density. Panel (a) uses 
the natural quantization axis where the skyrmion axes are 
aligned with the z direction, while (b) uses an arbitrary axis. 



The states produced from this melting are highly non- 
trivial, with strongly-correlated structures reminiscent of 
the multilayer quantum hall effect |3(j • Exactly how the 
textures discussed here are connected with the correlated 
states are a matter of current research. 
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the pseudospin-1/2 case at JILA and in the spin-1 
and 2 case at MIT 14] . An illustrative density profile 
for the three components of texture (H) is shown in fig- 
ure [21 Note that the Hamiltonian © is invarient under 
a global rotation of all the spins. In principle, this means 
that an experiment may create spin textures with a ran- 
domly chosen global orientation. Figure IT21 (a) and (b) 
shows the components of the same spin texture, but with 
the spins uniformly rotated. The differences between the 
two figures demonstrates how sensitive the component 
densities are to global spin orientations. 

A second detection scheme makes use of the birefrin- 
gent properties of a spinor condensate. As shown in |29| . 
light, when detuned from the fundamental transition of 
the atoms by a frequency which is large compared to the 
hyperfine structure but small compared to the fine struc- 
ture, interacts with the spin textures in a simple man- 
ner. At these frequencies, the ferromagnetic regions are 
optically active, meaning that circularly polarized light 
aligned with the spins travels at a different speed than 
the opposite circular polarization. Consequently, the po- 
larization axis of linearly polarized light rotates when it 
passes through a ferromagnetic region. The angle of rota- 
tion is proportional to the projection of the ferromagnetic 
order along the light propegation direction. No such ro- 
tation occurs when the light passes through a nematic 
region. 

Using this effect, one can envisage an experimental 
setup where the sample is imaged with polarized light, 
but with a crossed polarizer in front of the detector. Only 
light which has its polarization axis rotated by passing 
through a ferromagnetic region will reach the detector. 
Thus, one could directly image the ferromagnetic cores 
of the 7r disclinations found in the antifcrromagnetic gas. 
C. Quantum Hall physics 

At even higher rotation speeds, quantum fluctuations 
are expected to melt the regular lattices of spin textures. 



APPENDIX A: LOCAL VORTICITY 

In this appendix we derive equation i|15fl . which re- 
lates the curl of the velocity to the spin order parameter. 
Introducing a scaled order parameter <p c — xjj c /y/n, the 
superfluid velocity is given by 

v. = 5?-(#W c -0 c V#), (Al) 

where repeated indices are summed over. The vorticity 
is thus given by 

V x v s = — (V# x (A2) 
im 

This expression is related to l|15fl by noting that 

QabVQbc X VQoa = ftMctaVfi X Wc (A3) 

x V0 a (A4) 
+<t>* a <t> b <t>*<t>cV<i>i x V0 a (A5) 
+cj>*J b (t>tci>aV4>c x V0 C *. (A6) 

The terms (|A3|) and HA4[) vanish on account of respec- 
tively being antisymmetric in the indices b, c and a, a 
Noting th at 4> *V(fia = V(0*0 a ) - </> a V0*, one sees that 
the term (|A5ll is antisymmetric in the indices a, b and 
therefore also vanishes. Finally, using <j)* a 4> a = 1, one 
finds 

Qab^Qbc x VQ ca = -V0* x V0 C , (A7) 

which combined with (|A2|) yields l|15fl . Note that this 
result is not dependent on the atoms being spin 1 , but is 
completely general. 
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